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Abstract 



We find the exact non-perturbative expression for a simple Wilson loop 
of arbitrary shape for U{N) and SU{N) Euclidean or Minkowskian two- 
dimensional Yang-Mills theory (YM2) regulated by the Wu-Mandelstam- 
Leibbrandt gauge prescription. The result differs from the standard pure 
exponential area-law of YM2, but still exhibits confinement as well as in- 
variance under area-preserving difFeomorphisms and generalized axial gauge 
transformations. We show that the large N limit is not a good approximation 
to the model at finite and conclude that Wu's N = 00 Bethe-Salpeter equa- 
tion for QCD2 should have no bound state solutions. The main significance 
of our results derives from the importance of the Wu-Mandelstam-Leibbrandt 
prescription in higher-dimensional perturbative gauge theory. 
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QCD2, the dynamical theory of quarks and gluons in two dimensions, has been a fasci- 
nating testing-ground for some of the important concepts of high energy physics. A profound 
analysis of the theory was given by 't Hooft in the limit N ^ oo corresponding to an 
infinite number of colors, 't Hooft 's study has culminated in the explicit demonstration of 
quark confinement and in the numerical determination of the bound state meson spectrum. 

Technically, 't Hooft has worked in a generalized axial gauge, which means that the gauge 
potential is set to zero along a fixed direction. This direction is most conveniently chosen 
as light-like. Within this so-called light-cone gauge {A_ = 0), one of the two quark spinor 
components decouples, 't Hooft has treated the remaining infra-red singularity of the gauge 
field propagator in a way equivalent (cf B]) to a principal part prescription: 
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The highly consistent 't Hooft theory has been extended in many studies. In particular, 
the invariance of the meson spectrum with respect to a different choice of gauge, the space- 
like axial gauge, has been checked explicitly by a combination of analytical and numerical 
calculations 0. In this more complicated gauge one has to deal with both quark spinor 
components, and the integral equation for the meson bound states becomes two-dimensional. 
Nevertheless, an identical bound state meson spectrum was obtained [||. 

QCD2 would appear to be completely satisfactory were it not for two potential prob- 
lems, which turn out to be closely related. The first is, that using the principal part pre- 
scription Wick rotation to Euclidean space is impossible 0]. This was remarked early 
on by Wu P, who suggested to formulate QCD2 in Euclidean space. He considered the 
Wick-rotated kernel {k_ = ki — kQ —>■ ki + ik2) 
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and proposed to use symmetric integration. This means that the integrals with the kernel 
eq. (0) should exclude a small spherical hole of radius e around the origin of the (^1,^2) 
plane. 

Secondly, the principal part prescription for generalized axial gauges appears to he ir- 
reparably inconsistent above two dimensions since manifest perturbative renormalizability is 
lost. For light-cone gauges above two dimension, the regularization suggested by Mandel- 
stam and by Leibbrandt ^ appears to solve this problem (c/also 0). In two dimensions, 
this prescription gives 



1 1 V 



k"^ {k_ — it sgnA;_(_)^ k"^ 



— in S'{k-) sgn/c+. (3) 



Bassetto et al. noted that the Wu kernel (0) and the Mandelstam-Leibbrandt kernel 
are related to each other by a Wick rotation. 

Wu derived an effective integral equation for the bound state mesons in QCD2 at 
N = 00 with the regularization of eqs. (|^), (0). Physical observables should not depend on 
the regularization scheme, and it was widely suspected that the Wu-Mandelstam-Leibbrandt 
(WML) kernels should lead to the same results as the 't Hooft prescription. Past efforts to 
solve the Wu equation have however failed. 
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The apparent impossibility to resolve the meson bound state equation in the Wu for- 
malism is not a true impediment to study QCD2 in the WML prescription. Complementary 
insight can be obtained by probing the Yang-Mills field by static fiux lines instead of dy- 
namical quarks: the Wilson loops. This has already been the basic motivation of work by 
Bassetto et al. [Q. 

The present paper consists in a study of Wilson loops for QCD2. We first investigate 
perturbatively various contours in Euclidean and Minkowski space, up to 0{g^), using nu- 
merical and analytic calculations. We then derive the exact non-perturbative result. Our 
study confirms explicitly the finding that the WML regularization gives a result different 
from 't Hooft's. However, contrary to what was claimed before in the literature 0, the result 
fulfills the same consistency conditions as 't Hooft's Wilson loop: (i) it does not depend on 
the shape of the closed contour, but only on its area, (ii) the leading behavior at finite N 
is an area law behavior, indicating confinement, and (in) the result is gauge invariant in a 
(restricted) way which will be detailed later on. In our opinion, these tests put the WML 
regularization on an equal footing with the 't Hooft prescription. 

Before endeavoring on the detailed calculations in the WML regularization, we briefiy 
present the precise definitions and review the main results of the analogous calculation using 
't Hooft's prescription. Wilson loops are defined as 

Wc = (^PTr exp {tg I dx ■ A)). (4) 
iV Jc 

Here P denotes path ordering of the gauge field A along the contour C. In two Euclidean 
dimensions the standard, exact result can be obtained in a gauge invariant way for any num- 
ber of colors and by a variety of methods; e.g. using a manifest gauge invariant formulation 
like lattice gauge theory. For the case of a simple (i. e. not self-intersecting) U (N) or SU{N) 
Wilson loop Wc along a contour C one obtains 

where Ac is the area enclosed by C. Apart from the trivial factor N, the U{N) result does not 
distinguish between the Abelian and the non-Abelian case. Furthermore, the form of eq. 
agrees with the idea of a linear confining potential between sources. Wc depends on the 
contour C solely through the dimensionless combination g'^Ac- The general reasons for this 



fact were first emphasized in |Ty]: YM2 is invariant under area preserving diffeomorphisms 
rT| . Even in Minkowski space, the 't Hooft prescription (|I]) (or any other generalized axial 
gauge with principal part prescription) reproduces the exponential law (^ with Ac iAc 
where Ac is the Minkowski "area" enclosed by the contour. 

Let us first work in Euclidean space and use the prescription (Q). The path-ordered 
exponential in eq. (0) is defined as 



1 „ /"l rS2n-l 

^C = ^T.(-9T ds,x4s,)... dS2nX-is2n)^T{A+ix{s^))...A+{xis2n))), (6) 

N ^0 Jo Jo 



where x{s), s G [0, 1] parametrizes the closed contour C. The expectation value in eq. is 
to be evaluated by the Wick rule, and the basic correlator, obtained from eq. (@) by simple 
Fourier transform to configuration space, is {x± = xi =F 1x2) 
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= ^ klS,,, (7) 

Here we have also written out the matrix indices of the gauge field A^. 

It is simplest to consider a circular contour. Then the weighted basic correlator is inde- 
pendent of the variables s, s': X-{s) X-{s') — a;+(s'))/(a;_(s) — = (27rr)^. Now 
the integration over the path parameters Si, . . . , S2n becomes trivial, and the computation 
can be very simply performed, say up to three loops. One finds for U{N) 

= l-\ Ng'Ac + '/-N^ + I) g'Al - ^(^iV^ + l^iV) g^A^ + 0{g% (8) 

where Ac = vrr^ is the area of the circle. Evidently this three-loop result is incompatible 
with the standard law eq. (^) unless = 1. The technical reason is that the Abelian result 
is reproduced in the WML prescription by non-planar gluon exchange; as soon as > 1, 
crossing gluon lines acquire a special weight, as is evident from eq. (^: At 0{g'^), e. g., 
there are two planar and one crossed diagram, whereas at 0{g^) five planar and ten crossed 
diagrams are present. 

Is the result eq. (H), derived for a circular contour, generally true for all simple loops of 
the same area? Does eq. (||) remain valid for contours in Minkowski space, if we perform 
the analytic continuation Ac iAc"^ The answer to both questions is yes. The argument 



for the first affirmation is the invariance under area preserving diffeomorphisms |TT[], which 
should not depend on the regularization scheme. The second statement should be true since 
the kernels (|]) and (H) are related by analytic continuation, and since the Wu kernel (0) 
possesses a commutativity property already emphasized in the original work ||5|. 

Since the two above assertions are crucial for the following, we have convinced ourselves of 
their validity for a variety of simple contours such as ellipses, triangles, rectangles of various 
orientations both in Euclidean and in Minkowski space. Careful numerical evaluation of 
eq. (I) up to 0{g^) and in some case up to 0{g^) yielded in all cases perfect agreement with 
eq. to a level of precision of 10~^ . . . 10~^. Our numerical computations thus lend strong 
support to the formal general arguments mentioned above. 

Our findings are seemingly in disagreement with the analytical two-loop (i.e. 0{g'^)) 
calculations of Bassetto et al. for rectangular contours in Minkowski space [O- They 



consider two orientations of the loop. In the first case, the rectangle is oriented along the 



light-cone space and light-cone time axes. Agreement with eq. (g) can be established [13 
In the second case, the rectangle is oriented along the space and time axes. Here they 
concluded that the WML regularization violated the area law. However, the authors failed 
to realize that the complicated dependence of their final result on the aspect ratio of their 
rectangle exactly cancels, as is easily verified. If this cancellation is taken into account, the 
Bassetto et al. result is in perfect agreement with ours. 

We also investigated the invariance of Wc under transformations to generalized axial 
gauges (cf 0) 

(9) 



{ki+ik2Y {ki cos 6 + ik2 sin Oy 

for < 9 < 7t/4. In that case, symmetric integration amounts to cutting out a spherical 
hole in the {ki,k2) plane with ki = A;icos6' and ^2 = k2sm6. The propagator becomes 
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x+/x- with x± = X\l cos 6 =p 1x2/ sin 6. As before, the gauge invariance was 

checked numerically to a relative precision of about 10~^ for several values of 6. 

Having established that the choice of contour is arbitrary, we may restrict ourselves 
to the especially simple Euclidean circular contour. There we can find all further terms 
in eq. (|^), since the integrand in eq. (|^) is constant. The problem of determining the 
Wilson loop reduces to the purely combinatorial problem of finding the group-theoretic 
factors corresponding to the Wick contractions. Fortunately, these factors are generated by 
a simple matrix integral: 

Wc = ^ JVF exp{- ^TtF^) ^Tr exp (tg^F). (10) 

For U{N), VF denotes the flat integration measure on the space of hermitian N x N 
matrices: VF = Uf=idFiiUf<j d{ReFij)d(lmFij). Z is a normalization factor: Z = 
JVF exp ( — ^TrF^). This matrix integral has been evaluated with a variety of meth- 



ods [Tj]. The final result constitutes the exact expression for the Wilson loop at any A^: 



,u{N) _ ( I 2 . . ^ f dz , ^ , Jz^\\^ 



y^c' =^^v{--/Ac)-j—^xv{-9'Acz)[-^) . (11) 

The contour integral, which encloses the multiple pole at 2; = 0, gives a Laguerre polynomial 
in g^Ac of order — 1: L\q _^{^g^ Ac) . The first few examples are: 

We^W = e-^^^-^^ = (1 - i,Mc)e-^^^^S W.^^^^ = (l - /A + 

(12) 

For SU{N), the integration measure in eq. ([I0|) has to be modified to enforce tracelessness: 
VF —>■ VF 6{TtF). The U{1) part decouples from the U{N) Wilson loop and we find 

Wr(-)=exp(^/A)Wr^ (13) 

Our result coincides with the usual expression (||) only in the Abelian U (1) case. Incidentally, 
this explains why Wu and Stamatescu were able to reproduce the standard solution of 
the Schwinger model (i.e. two-dimensional QED) with the WML regularization. 

We now study the exact formulas ([Hp, (|T3p in 't Hooft's large A^ limit: N ^ 00, g ^ 0, 
with fixed g^ = Ng"^ . Eqs. (|ll|)! (HH) become the integral representation of a Bessel function: 



W^^^^ = W^""^^^ = Ji{2~gJAc). (14) 



Curiously, the factors exp(— |(7^^c) present at all finite A^ (cf eqs. (pA]), (|12D) disappear in 
the N ^ 00 limit. But it is precisely these factors which lead to confinement at finite A^. 
Very much unlike the usual theory, where from eq. (^ we have W^^""^ = exp{—^g'^Ac), in 
the present model the behavior changes qualitatively at A^ = cxd, as Ac —>■ 00 

W^^°^^ ~ -^ifAc)-'^ cos(2^V^ - (15) 
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This fall-off at large areas Ac is far too slow to ensure confinement! We therefore conclude 
that Wu's model of mesons has to be studied at finite N; the Bethe-Salpeter equation written 
down in is not expected to lead to a discrete meson spectrum. 

What have we achieved so far? We have presented evidence that, besides the 't Hooft 
principal value prescription, the WML symmetric integration yields a consistent yet different 
theory. We have found this theory to be qualitatively different in the ^ oo limit. This 
allows us to solve the old riddle of why bound states for the Wu equation have never been 
found. It also allows us to expose a case in which the N ^ oo limit is invalid. 

At finite N, we are lead to conclude that two physically different, consistent, gauge 
theory formulations in D = 2 exist, unless inconsistencies of the WML prescription end up 
being found on a subtler level. It is indeed known that QCD2 is special in that it can be 
considerably generalized by adding higher powers (products of traces) of the field strength to 
the action; such terms are irrelevant in Z) > 2 but begin to scale at exactly D = 2 |10[, [ll6| . 
Therefore, as already pointed out in |]TB[, continuing a D-dimensional formulation of QCD 
down to D = 2 might not give the usual TrF^ theory, but one of the generalized theories 



(gQCD2) |T^. However, this does not explain why there are coexisting theories with the 
same Lagrangian C = -\TtF^. It is easily seen that the WML prescription in Euclidean 
space cannot be described by a standard (even generalized) lattice gauge theory. It would 
be very interesting to find a lattice discretization for this theory. 

We emphasize that the two prescriptions naturally result from different points of view: 
the standard prescription is obtained from a number of strictly two-dimensional quantization 
methods (Euclidean lattice gauge theory, light-cone quantization) while the WML prescrip- 
tion seems compelling from the viewpoint of higher dimensional Minkowski perturbative 
QCD. It appears that "theory space" is big enough in D = 2 to allow coexistence of these 
theories. 
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